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Abstract 



We derive the exact semiclassical expression for the second inertial parameter B for the 
superfluid and normal phases. Interpolation between these limiting values shows that the 
function B{I) changes sign at the spin Ic, which is critical for a rotational spectrum. The 
quantity B turns out to be a sensitive measure of the change in static pairing correlations. 
The superfluid-to-normal transition reveals itself in the specific variation of the ratio B/A 
versus spin / with the plateau characteristic of the normal phase. We find this dependence 
to be universal for normal deformed and superdeformed bands. The long plateau with a 
small value B/A^ A~^^^ explains the extreme regularity of superdeformed bands. 
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1 Introduction 



Recently phase transitions in mesoscopic systems have been a subject of intense discussions 
in nuclear and solid-state physics. Isai Isidorovitch Gurevich appears to be the first to bring 
up (in 1939) the concept of the temperature phase transition to nuclear studies His 
prediction was based on the observation that the level density of the resonant states formed 
by a thermal neutron capture is an unsteady function of the atomic mass number with a 
maximum in the rare- earth nuclei. 

The problem of the rot at ion- induced transition from the superfluid to the normal phase 
in nuclei has been a foremost theme in high-spin spectroscopy since Mottelson and Valatin 

predicted a pairing collapse in rapidly rotating nuclei. This effect can be understood by 
an analogy with a superconductor in a magnetic field. In a deformed nucleus the Cooper 
pair is formed by two nucleons with opposite single-particle angular momentum projections 
±m. Being time noninvariant (as a magnetic field) the Coriolis force in a rotating nucleus 
acts on both nucleons in opposite directions and tries to decrease the spatial overlap of 
these time reversal orbits. The Coriolis force increases proportionally with the spin of a 
band. Therefore, at some critical spin one may expect that all pairs are broken and pairing 
correlations disappear completely. The phenomenon can be observed by the crossing of the 
ground state superfluid band with the band based on the normal state. Thus, the rigid body 
moment of inertia corresponding to the second band appears to be the obvious signature of 
the pairing phase transition. 

However, this regime of the transition to the normal phase is not realized in nuclei because 
they are finite systems with a shell structure and a small number of nucleons involved in 
pairing correlations. The Coriolis force in a rotating nucleus is proportional to the value 
of the single-particle angular momentum j of a nucleon. Thus, the Coriolis antipairing 
effect is strongest for nucleons occupying the states of subshells with the largest j. In the 
vicinity of the Fermi surface, these so-called intruder orbitals arise from the j = M + 1/2 
subshell, where M is the principal quantum number of the shell above. Therefore they are 
distinguished from other states of the unfilled shell by the parity. At normal deformations 
(ND), the intruder states retain their j quantum number, while at super deformations (SD) 
the j-subshell notation becomes less appropriate due to mixing. Initially the Coriolis force 
breaks only the one Cooper pair that belongs to the intruder orbitals, whereas the rest of the 
pairs stay correlated. The band built on such two-quasiparticle excitation (the rotationally- 
aligned band) is characterized, due to the blocking effect, by appreciably smaller pairing 
correlations than in the ground state band. Having the largest moment of inertia, the former 
crosses the ground state band and becomes yrast. The relevant phenomenon is backbending. 
Subsequent breaking of correlated pairs and their alignment makes the internal structure 
of the yrast band non-homogeneous and the transition to the normal phase configuration 
dependent. 

The standard definition of the phase transition is based on the mean-field approxima- 
tion in which different phases are distinguished by the order parameter, i. e., the static 
pairing gap A. However the mean field approach to the nuclear pairing correlations faces 
a fundamental problem of the quantum fluctuations, which become quite strong for finite 
systems. The fluctuating part (dynamical pairing correlations) of the order parameter 
is comparable with A in the transition region. The fluctuations smear out a sharp phase 
transition and make it a difficult issue to find the experimental signature of the superfluid to 
normal phase transition in rotational bands. For example, the dependence of the kinematic 
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dynamical Q"*^^-' moments of inertia on the spin / is not a definite indicator of the 
phase transition. Experimental evidence of the pairing phase transitions has been discussed 
usually in terms of the relative excitation spectrum. As shown in Ref. the disappearance 
of static pairing leads to a change of the excitation spectrum, from the quasiparticle to the 
particle-hole spectrum. Unfortunately, the application of this criterion to normal deformed 
(ND) bands |]^, ^, ^ shows that this method is not free from ambiguities. 

Meanwhile, it is well known that the change in the system internal structure manifests 
itself in the modification of its collective excitation. The examples for finite quantum systems 
are numerous. A classical one is the transition from deformed to spherical nuclei. In this 
case, the rotational-vibrational spectrum transforms into a pure vibrational one. The study 
of the bifurcations in rotational spectra [^] shows an intimate connection between internal 
and rotational motion. For example, the angular momentum alignment in a band (the change 
of a coupling scheme) is observed as an increase of the energy signature splitting 

The transition we study is more delicate. Consider the simplest rotational sequence with 
the parity and the signature (vra) = (+0). The relevant energy spectrum can be parametrized 
as follows: 



E{I) =AI{I + l) + Bl\l + lf 



where A = h'^/2Q^^^ and B are the first two inertial parameters. The spectrum ( p..l|) under- 
goes a noticeable modification if, for example, the second inertial parameter changes sign at 
some spin Ic- For I < 1^ the spectrum is compressed relative to the rigid rotor spectrum 
because B is negative for the low-/ states. However, for I > Ic the spectrum becomes ex- 
tended. The effect can be visualized by using the / dependence of the ratio B/A. The main 
objective of the present paper is to study this dependence. 

The parameters A and B are determined by the 7-ray transition energies E^[I) = E{I + 
2) - E{I) as follows: 
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The coefficient B characterizes the nonadiabatic properties of a band and is very sensitive 
to its internal structure. It also realizes the relationship between kinematic and dynamic 
moments of inertia. Using the well-known expressions for these quantities (see, for example, 
Ref. IQ]) and the last formula ( |1.2| ) we get 
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Thus, the parameter B{I) is proportional to the difference '^^^^ — ^S^^^ in the high-J limit. 
The ratio B/A also determines the convergence radius of the rotational energy expansion in 
terms of /(/ -|- 1) 0. Faster convergence is obtained with the Harris formula 



E(uj) 



;i-4) 



which is based on the fourth-order cranking expansion 
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where 

p{n) is 

the nth correction to the nucleus density matrix; is the projection of the 
single-particle angular momentum operator onto the rotational axis x, which is perpendicular 
to the symmetry axis z\ and uo is the rotational frequency. The latter depends on the system 
angular momentum and is determined by 



n^I{I+l) =auj + (3uj^ + ... . (1.6) 
It follows from Eqs. (0), and (|OD that 

The problem of the microscopic calculation of the parameter B for ND nuclei has attracted 
considerable attention (see the review article |]10| and references therein). It has been shown 



that this quantity receives the contributions from four types of nonadiabatic effects: 

(i) perturbation of quasiparticle motion by rotation (quasiparticle alignment), 

(ii) attenuation of pairing correlations by the Coriolis force (Coriolis antipairing effect), 

(iii) a change in the deformation of nuclear self-consistent field (centrifugal stretching 
effect), 

(iv) vibration-rotation interaction. 

The first attempt to estimate B was made by the author together with Grin' [11|. A Green's 
function formulation of the Hartree-Fock-Bogolubov (HFB) method was used to find the 
expansion ( |1.4| ) for the axially deformed oscillator potential as the self-consistent field. It 
was shown that the first and the second effects yield B/A ~ A~^^^ while the centrifugal 
stretching contribution is A^^^ times smaller for well deformed nuclei.^] In the subsequent 
work [1^, the author found that the vibration-rotation contribution to the parameter B 



accounts for the same A^''^ fraction of the main effects. These results were confirmed by the 



calculations of Marshalek |T3| with the more realistic Nilsson potential. 

Thus, the first two effects are dominant for well deformed nuclei. The quasiparticle 
alignment depends strongly on pairing correlations because the pairing force tries to bind 
pairs of particles in time reversal states, reducing the ability of nucleons to carry an angular 
momentum. Therefore, the parameter B is very sensitive to the variation of the pairing 
correlations along a band. 

One of the amazing features of SD bands is the extreme regularity of their rotational 
spectra. To demonstrate this feature, the rotational spectra of different axial systems are 
compared in Fig. 1 with their rigid rotor counterparts. The comparison shows that the SD 
band ^^^Pb(l) is more regular than the ND band of ^ss^j g^^^ even the band of the simplest 
H2 molecule Having the ratio B/A ~ 10~^, the band ^^^Pb(l) is not a champion among 
SD bands. For ^^^Dy(l), the ratio is of the order 10"^ and this is 1000 times smaller than 
the above estimation B/A^ A~^^^. Thus, an SD nucleus is the best quantum rotor known 
in nature. Although numerous theoretical calculations (see, e.g., |15|, |1^, [1^, |I^) successfully 



reproduce the measured intraband 7-ray energies, the underlying microscopic mechanism of 
this phenomenon has yet to be well understood. 

In this paper we will reveal an interconnection between the extreme regularity and the 
transition from the superfluid to the normal phase. The key to our theoretical approach 
lies in the calculation of the second inertial parameter. Compared to ND bands, there are 
two features of the pairing correlations in SD ones which prevent us from using the results 

""^We use the estimation A ^ £f^~^^'^, where ef is the Fermi energy and A is the mass number. 
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of previous theoretical calculations of the parameter B for superdeformation. First, due to 
the large shell gap stabilizing the SD minimum, the static pairing field A is small and can 
be commensurate with its fluctuation 6 A. Qualitative conclusion concerning the role of the 
static and dynamic pairing in SD bands is presented in Ref. |[19[. Second, since intruder 
single-particle states, which are unavailable at normal deformations, appear near the Fermi 
surface in the case of superdeformations, it is necessary to go beyond the commonly used 
monopole pairing force pO[. The gauge invariant pairing interaction expands the correlation 
space and stabilizes the pairing field. The coordinate dependent (nonuniform) pairing is also 
crucial for conservation of a nucleon current in a rotating nucleus ||21|| . 

To avoid calculation of the parameter B in the transition region, where pairing fluctua- 
tions play an important role, an interpolation between the values Bg and Bn is used. The 
former is associated with the superfiuid phase (where A 3> 6 A) and the latter is related to 
the normal one (A = 0). Thus, pairing fluctuations are unessential for these regions and we 
can use the mean-field approach. In the calculation of Bg, the nonuniform pairing induced 
by rotation is taken into account by using the method of Ref. |TT[ . It should be noted that 
the quantity Bs found in the cited work is inapplicable for superdeformation. 

The paper is organized as follows: In Section ^ the basic equations of the cranked HFB 
theory are presented in the framework of the Green's function method. The spinor form of 
the Gor'kov equations is used to simplify calculations in the higher orders of the perturbation 
theory. In Section |^, the exact expression for the second inertial parameter in the superfiuid 
phase is derived with this technique by applying a semiclassical approximation. The result is 
valid for an arbitrary nuclear mean field. The exact analytical expression for Bg is obtained 
in Section ^ in an axially deformed oscillator potential. In this section we also consider some 
limiting cases for this quantity. Of special interest is the limit of noninteracting nucleons. It 
is shown that the relevant parameter i3„ is positive and smaller than Bs- The comparison 
with available experimental data for SD and ND bands is presented in Section ^ Section |^ 
concludes and summarizes the paper. The preliminary results of the present work have been 
published in Refs. P^, 



2 Green's function formalism in the cranking Hartree- 
Fock-Bogolubov method 

2.1 Cranked Gor'kov equations 

Our consideration is based on the shell-model Hamiltonian consisting of the cranked one 
body term 

h^{r) = -^ + U{r)-uj-i, uj{uj,0,0} (2.1) 

(where p and M are the impulse and the mass of a nucleon respectively), and the residual 
short-range interaction, which is specified by the two body delta-interaction 

v{r,r') = -g6{r-r'), g > 0. (2.2) 

For simplicity, we neglect the spin in the cranking term and consider only the orbital part £ 
of the angular momentum j. We will also neglect a weak dependence of the self-consistent 
deformed potential U on rotation (centrifugal stretching effect). 
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In the coordinate representation, the Gor'kov equations have the form 
[e - K{v) + £^]G(r, r', e) + A(r)F+(r, r', e) = 5{v - r'), 
[e + /i+(r) - £i.]F+(r, r', e) + A+(r)G(r, r', e) = 0, 

A*ir) = gi^F^{v,v,e). (2.3) 

c 

The functions G{r,r',e) and F"''(r, r', e) are the Fourier transforms of the Green's functions 

G(r, r', t-t') = -2(<l>^|r{^(r, t)^+(r', t)}|<l>^), 

F+(r,r',t-t') = -z($^+2|T{^+(r,t)^+(r',t)}|<|.^)e-2-^*, (2.4) 

where $Ar and are the eigenfunctions of the ground state of a system of and + 2 

interacting particles, ?/'"'" and ip are creation and annihilation operators in the Heisenberg 
representation, T is the time ordering operator, and Ef is the Fermi energy (the chemical 
potential of the system). The contour C consists of the real axis and an infinite semicircle 
in the upper half-plane. 



2.2 Properties 



In obtaining Eqs. (|2.3| ) the particle number nonconserving approximation has been used. 
In a spirit of the mean-field approach we neglect a difference between the functions and 
Thus, the Gor'kov equations describe a system with the broken gauge symmetry 
associated with the particle number. However, the average particle number is fixed. This 
is achieved by adding the term —epN to the Hamiltonian. The Lagrange multiplier ep is 
determined by the equation 

N = [ dr <f ^G{r,r,e). (2.5) 



2m 

c 



Equations (|2.3|) are also noninvariant with respect to the more general gauge transfor- 



mation (the local Gallileian transformation |^ 



V'(r,t)-^^(r,t)e^^W, (2.6) 

where 0(r) is an arbitrary function of the space coordinates. The quickest way to show this 
is to introduce the vector potential A = [co'r]/2 that allows us to rewrite the Coriolis term 
V = —uj ■ i in the form — 2p ■ A. It is seen that the Hamiltonian (^TTJ) lacks the term 2MA^ 
which is absolutely necessary for the gauge invariance of Eqs. ( p.3|) . However, since the 
two-body interaction (|2.2|) is invariant under the Gallileian transformation the conservation 
of nucleon current is ensured. 

The current density is expressed in terms of the Green's function G as follows 



c 

With this definition, we find 



divj(r) = hm f T^|A„(r) - A:(r')]G(r, r', e). (2.8) 

C 
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Using the first of Eq. ( p^.3| ) and their complex conjugate equation 

[e - ht{v') + eF]G{v, r', e) + A*(r')F(r, r', e) = 5{v - r'), (2.9) 

we finally obtain 

div j(r) = j ^[A(r)F+(r, r, s) - A*(r)F(r, r, e)]. (2.10) 

c 

The right-hand side of this equality vanishes due to the third Gor'kov equation. Because the 
latter is derived assuming delta-interaction we should conclude that the form of a two-body 
interaction is essential for obtaining the current conservation. In particular, the commonly 
used monopole pairing interaction is not invariant under the transformation ( |2.6|) . Therefore 
it does not conserve the current density in a rotating nucleus. The case of an arbitrary pairing 
interaction is considered in Ref. [p5|] . 

2.3 Matrix form of the Gor'kov equations 

The two-dimensional form of the Gor'kov equations is very convenient to use in our calcula- 
tions. Let us introduce the second pair of the Green's functions G^{t, r', e) and -F(r, r', e). It 
is easily proved |]TT| that the four equations for these functions can be written in the matrix 
form 

( A* e + h+-eF\( F G ^ 
\e-K + eF A )\G+ F+ ^ 

The operator involves the real h and the imaginary V parts. The former is the Hamilto- 
nian of the deformed mean field, whereas the latter denotes the cranking term, V = —Loix- 
Separating the quantity A into real and imaginary parts, 

A = A + A, A* = A-A, (2.12) 

we can rewrite Eq. ( ^.llj ) in the compact form: 

[^p + A- a,V - asA)G{r, r', p) = 6{r - r'), (2.13) 

where G is the matrix of the functions G and F involved in Eq. ( p.ll| ), p = aip + a2{h — eF), 



i(5(r-r'). (2.11) 



(Jq are the Pauli matrices, and p = —ie. We omit the unit matrices before the terms with A 
and 6{r — r'). The functions G and F can be written as traces of G in the following way: f\ 

G(r,r',e) = ^Sp{(ai-z<T2)G(r,r',p)}, 

F(r, r', e) = ^Sp{(l + as)G{r, r',p)}. (2.14) 
Therefore, the equation for A(r) is 

H^)=gj ^Sp{(l + a3)G(r,r,p)}, (2.15) 

c 

and the one-particle density matrix of the system is given by the expression 

p(r) = j ^Sp{(ai-2a2)G(r,r,p)}, (2.16) 

c 

where the contour C is obtained from G by the 90° rotation. 

^We use the symbol tr in the space of single-particle states of the Haniiltonian h, the symbol Sp in the 
spinor space, and Tr in the combined space. 
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2.4 Perturbation theory 



We apply the Green's function method to calculate the parameter B. As follows from Eqs. 
( |1.5|) and (|1.7| ), this requires the perturbation theory of third order in the cranking term V. 
According to Ref. [0, a considerable computational effort is needed to mold the result into 
a tractable form. The matrix representation of the Gor'kov equations allows us to elaborate 
on the elegant form of the perturbation theory which considerably simplifies the calculations. 

We now proceed to treat Eq. ( |2.13| ) by the method of successive approximation. Con- 
sidering V as a weak perturbation, we expand the Green's function and the self-consistent 
quantities A and in the powers of a small parameter: 



G — Gq + Gi + G2 + Gs + 



A = A(°) + A(i) + A(2) + A(=^) + 6^ = eP + ef + ... . 



-(2) 



(2.17) 



The nth order corrections to the last two quantities are determined by Eqs. (|2.15|) and (|2.5| ), 
respectively. Since V = —V*, the corrections of odd order to A are purely imaginary and 
those of even order are real. The effect of the second order correction to sp gives a negligible 



small contribution in the second inertial parameter |T3[. Thus, we will use the zero-order 
approximation for this quantity. 



It is natural to work in the basis of the eigenfunctions of the Hamiltonian without 
the cranking term, 

{h - eF)Vu{r) = Pu^A^), (2-18) 

where Pi, is the energy e^, of the single-particle state u relative to the Fermi energy, p^, = 
Ey — Bp- In this basis, the rath correction to the Green's function has the form 



G„(r,r',p) = 5:GS(p)^.(r)^:,(r' 



(2.19) 



First of all, we find the solution of the unperturbed equation ( p.l3|) , 

(^p + A)Go(r,r',p)=5(r-r'), (2.20) 
with the constant pairing gap A^^^ = A. Substituting Eq. (|2.19| ) into (|2.20|) , one finds 

.(o)/_ N iPv-^ 



Gyly'iPu) 



(2.21) 



where Pu{p,Pi^) is the two-dimensional vector and p^, = crip + d'2Pu. The gap equation ( ^.151 ) 
takes the simple form 



2E, 



pI 



A2 



(2.22) 



The equation has the solution A = const for the self-consistent potential with a flat bottom. 
In obtaining Eq. ( |2.22| ) as well as in subsequent calculations, it is essential to compute the 
traces of the products of the Pauli matrices. We can readily see that the trace of an odd 
number of matrices o"i and o"2 vanishes and that of an even number is given by the expressions 



-Sp((7a(T_a) — 5ap, -Sp((TaO-/}a^(T5) — 6ap5^s — Ea/S^-rS, 



(2.23) 
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where e = ia2 is a fully antisymmetric matrix. 

In the first order, Eq. ( p.l3|) involves the two perturbing terms V and A*^^); 

A)Gi(r,r',p) = [ail^ + a3A«]Go(r, r',p). (2.24) 

The solution of this equation is obvious: 

Gi(r,r',p) = J Goir,q,p)Wiq)Go{q,r',p)dq, (2.25) 

where the operator 

W = aiV + agA^^) (2.26) 

is introduced to write the corrections to the unperturbed Green's function in the simple 
symbolic form: 

Gi = GqWGo, G2 = GoWGoWGo - GqA^^^Go, 

G3 = GoWGoWGoWGo - GoWGoA^^^Go - GoA^^^GoWGo + Go'XsA^^)^^. (2.27) 

Here the integration over intermediate coordinates q is implied. Using these formulas, one 
can prove by straightforward calculations the following identities: 



I ^Sp{G2m(r,r',p)} = 0, J ^Sp{aiG2.(r, r',p)} = 0, 

c c 

dp f ^ ^ , [dp 



I /Sp{a2G2m(r,r',j9)} = 0, / /Sp{a3G2.(r, r',p)} = 0. (2.28) 

c c 

In order to find the self-consistent solution, we have to show how A*^""-* is obtained from 
Gn- We will consider the unperturbed pairing gap A as a parameter of the theory. This 
allows us to eliminate the interaction constant g. Multiplying the zero-order equation ( p.22|) 
by A'^"^(r) we write the result in the symmetric form: 

AA^-\v) = gf ^Sp[Go(r,r,p),A(")(r)] + , (2.29) 
J on 
c 

where [...]+ is the anticommutator of corresponding operators. With this ansatz, the integral 
equation for even order corrections is given by: 

/ ^Sp{2AG2^(r,r,p) - [Go(r, r,p), A(2^)(r)] + } = 0, (2.30) 
J Ztx 

c 

while that for odd order ones has the form: 

/ ^Sp{2Aa3G2m(r,r,p) - [Go(r, r,p), A(2^+i)(r)] + } = 0. (2.31) 

c 

Similarly, eliminating g from the equations for A^") and A*^"^^) (n is odd) yields 

J ^Sp{A(")(r)G„_i(r,r,p) - A("-^)(r)a3G„(r, r,p)} = 0. (2.32) 

c 



In the same manner, with additional integration over dr, we can derive 



f -^Tr{A(")a3G„_2(p) - A^'"''>asG,,ip)} = 0. (2.33) 

J ZTT 

C 

Here Tr refers to the spinor and to the single-particle spaces simultaneously. 

The solution of the HFB equations presented in the form or the successive approximations 
by the formulas ( |2.21| ), ( |2.22| ), ( p.27| ), ( |2.3CI| ), ( p.31| ) takes into consideration nonuniform 



pairing correlations induced by rotation. The nonuniform pairing originates in higher orders 
of the perturbation theory, while the nonrotating system is approximated by the constant 
pairing field. In spite of an apparent eclecticism, our solution does not violate the current 
density conservation, divj = 0. This is obvious in the zero-order approximation. In the 
even orders of the perturbation theory divj^^*-' vanishes due to the first and four identities 
of (|2.28| ). The odd corrections to the quantity ( |2.10[ ) can be transformed into the expression 



divj(2^+i)(r) = ^ / ^Sp{A(2'+i)(r)G2.-2Kr,r,p) - A^^O^gG^.-amlr, r, p)}, (2.34) 

which vanishes due to the equation ( p.32| ). Thus, we have found the consistent solution, 
which is certainly more general than initially supposed. It can be used for an arbitrary 
single-particle potential without danger of coming to contradiction (see also Ref. |2^). We 
will apply this solution to calculation of the second inertial parameter. 
Finally, the following identities are useful in calculations: 

2A / gTr{A(2^)G2.(p)} = / ^Tr{A(2^)[A(^^), Go(p)] + }, (2.35) 
c c 

2A / gTr{A(2^+i)a3G2.+i(p)} = / ^Tr {A^'^^'^ [A('^+'\ G,{p)U}. (2.36) 



c c 



These identities are obtained by multiplying Eq. ( p.30| ) and ( |2.31| ) by A*^^*) (r) and A*^^*+^^ (r), respectively, 
and integrating over dr. 



3 Calculation of the second inertial parameter 
3.1 General expression 

As follows from Eq. ( |1.7|) , the derivation of the parameter B is reduced to the calculation 
of (3. This latter quantity is convenient to deal with. The third-order correction to the 
density matrix can be obtained from Eqs. (|2.16| ) and ( p.27|) . If we take into account the 
third equation from ( p.28| ), substitute 

ai4 = --(W^-'^3A«), (3.1) 

UJ 

and use the identity ( |2.33| ) with ri = 3, we get after some simple algebraic calculations 

A = ^tr{4p'»'} = -ij/ |Tr{H/G3(p)} + ^/ !^Tr{A':%,GM}. (3.2) 

c c 
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where the subscript s means that the relevant quantity refers to the superfluid state. Here 
and below we will use A*^^*"*"^) {i = 0, 1) instead of A'^^*"''-'^^ Fortunately, the terms with A^^'' 
are eliminated from ( p.2|) after inserting the expression for G3 (|2.27| ). The resulting formula 



involves the corrections A*^^) and A*^^) only. It is convenient to transform the terms with A*^^-* 
into a quadratic form in this quantity. Referring to the definition of the function G2 from 
Eqs. ( |2.27| ), we find that the term involving A'-^-' becomes 

^Ti{A^'^ Go{p)WGo{p)WGo{p)} 
Zn 



c 



= J gTr{A(^)G2b)} - / ^Tr{A(2)Go(p)A(2)Gob)}. (3.3) 
c c 

With the help of the identity {^^) it is easy to show that the first term on the right also 
transforms into a bilinear form of A^^-*. Combining this result with ( |3.2| ) we obtain the final 
expression for the parameter Pg-. 

Ps = -—J ^Tt{WGo{p)WGo{p)WGo{p)WGo{p)} 



c 



+ £Tr{2AA(2)Go(p)A(^)Go(p) + A^'^Goip), A^^'l+j. (3.4) 

c 

This is an exact formula for the calculation of the second inertial parameter. The first term 
describes the joint effect of the Coriolis force and the nonuniform pairing field A*^^)(r) on 
independent quasiparticle motion. In the limit of the monopole pairing interaction, which 
corresponds to the uniform pairing field (A*^^-* = 0), this agrees with the term Be found by 
Marshalek [T^. The second term arises only due to the modification of the pairing. In the 



limit of uniform pairing this term describes the Corriolis antipairing effect (see Appendix 
A). 

3.2 Semiclassical approximation 

To proceed further we should find the first and second order corrections to the pairing field. 
As shown in Appendix A, the solutions of the corresponding integral equations are 

A«(r) = -^A4, A(^)(r) = ^D,e (3-5) 

where Di and D2 are the amplitudes, which are found in a self-consistent way. In order to 
learn more about the nonuniform pairing, we suppose that the self-consistent potential of a 
deformed nucleus is of the form 

UM^U{^^^^;). (3.6) 

where a and b are the half-axes of a nuclear spheroid. Then we obtain 

dll dll h'^ - 

^^ = y-^-^-^ = —i^y^u\p). (3.7) 
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Thus, the rotation induced pairing field is a function of the space coordinates only. In the 
first order in rotation, the nonuniform pairing field is proportional to the spherical harmonic 
y2±i- This was a motive for introducing the quadrupole pairing (see and references 



therein). The second correction A*^^^ shows, however, that higher multipoles are also involved 
in nonuniform pairing. 

Using the expression for A^^-* and the obvious formula for matrix elements M^^^, = i{py — 
Pv')^vu'i can represent each of 16 sums in the first term of Eq. (|3.4| ) in the general form 



\ ^ nx nx nx nx 

f dp Q4{p]Pi,P2,P3,Pa) „^ 

2.^12^34^41 / 2^ + ^2)(p2 + ^2)(p2 + ^2)(p2 + ^2) ' l^'^i 



C 



where the summation indices 1,2,3,4 refer to the single-particle states u with the energy pi, 
of the Schrodinger equation ( |2.18 ) and Q4 is a polynomial of the fourth power in {p,Pu) and 



Di, which is derived from calculating the trace of the product of the Pauli matrixes and 
factors {ipi, — A) . 

To evaluate this sum we use the method proposed by Migdal |]2T[. Let us note first that 



at the fixed state 1 the indexes 2, 3, and 4 take only few values permitted by the selection 
rules for the matrix elements of the operator i^. After integration over dp we obtain the 
function of the variables pi, ...,p4. When considered as the function of the variable pi, it has 
for fixed differences pi—p2,Pi—P3, and pi—pi a sharp maximum at the Fermi surface with 
an approximate width A. Since the average level spacing for ND nuclei, 5e ~ Sp/ A, is very 
small compared to A ~ ep/ ^^^^ a large number (~ A^^^) of levels fall within the interval 
A. For this reason we can make with semiclassical accuracy the following substitution in the 
sum (pD: 

dp_ Qi{,p;Pi,P2,P3,Pi) 

27r (p? + A2)(p2 + A2)(pi + A2)(p1 + A2) 

r/ _ X [ <^ Q4{P;PUP2,P3,P4) .„qx 

^ ^ J 2n (p? + A2)(pi + A2)(pi + A2)(p1 + A2) ' ^^'^^ 
Similarly, the second term of Eq. (|3.4D can be approximated by the expression 



^ ^Ia(2)|2 f dpi {pi-p2f+^^ 



47r (p2 + A2)(p2 + A2 

When calculating these integrals over dpi = dpdpi, it is convenient to use the Feynman 



covariant integration method because good convergence of the integrals allows the inte- 
gration over dpi to be extended from —00 to 00. For details of calculations see Appendix 
B. The integral in ( |3.9| ) depends on the three independent differences pi — P2, Pi — P3, and 
Pi — Pa- To represent the final result of the semiclassical approximation in a symmetrical 
form, we introduce the six energy differences p,^,^' = Pu — Pu', ^ < i^' ■ Collecting all the 
integrals of the first and the second terms in ( p.4|) , we find 

/?s= ^^X!^^2^23^34^4i^(3^i2,a;23, 2:34, 2:41; a;i3,X24)5(£:i-£:F), (3.11) 

where 5 function denotes that summation over the state 1 is substituted, according to the 
semiclassical approximation, by integration over its quantum numbers. 
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The function F, which depends on the six dimensionless differences Xypi = {eu—£u')/'^'^, 
is divided into two parts: 

F = f {Xu, X23, XSA, Xi3, + SDlxi2X23X34.X4lh{xi3) . (3.12) 

The first one is relevant for the first term of Eq. (3^). It is convenient to represent this 
function in the form 

/= -(1 + A + A + p3)G+il + Pi)H, (3.13) 



where 



G = — ^^'^illl |(l-Dia;i2) l+Xi2+a:23a^4i--Di [0:23(1 -a;i2X24) 

+ X4;^(l+Xi2a;i3)J — -Dx'^233^342:4l(2:23 + 2;4l) ~ -D;^Xi2X23X34X4;^ 

- Di(x34 - -010:12X23X41) (X34 + X12X13X24 - -01X12X23X41)} (3.14) 



and 



H-- 



h{x 
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X12X23X34X41 
Here the functions 



1-L)l(x^2 + X23 + X34 + X4i) + L)i(xi2X4i+X23X34)^ . (3.15) 



g{x) = h{x) = (1 + x')g{x) (3.16) 

XV 1 + x^ 



are associated with the Migdal moment of inertia The amphtudes Di and D2 of the 
nonuniform pairing field are determined by Eqs. ( |A.4| ) and ( [A.8| ), respectively. The operator 
Pi permutes the indices u of single-particle states in all the dimensionless differences, on 
which the functions G and H depend. When applied to x^^', we get 

PiXu,u' — Xjjj^iyj^i^ (3.17) 

subject to z/mod4 = v. It is easy to prove the following symmetry properties of the function 
F: 

PiF = P2F = P3F = F, 

F{xi2, X41, X34, X23; -X24, -X13) = i^(Xl2, X23, X34, X41; Xi3, X24)- (3.18) 

The above formulas give the semiclassical expression for the second inertial parameter in 
the superfiuid phase. The solution takes into account the effect of rotation on the Cooper 
pairs in the gauge-invariant form. The result is expressed entirely in terms of matrix elements 
and corresponding energy differences providing the constant pairing gap A is fixed for a 
nonrotating nucleus. It is valid for an arbitrary nuclear mean field with a stable deformation. 
This allows one to study an interplay between rotation, pairing correlations, and mean field 
deformation in ND and SD bands. 

We first estimate the quantity jSg and find the small parameter of the perturbation theory 
we used. To get an estimate of Ps for ND bands, we observe that the matrix element i^^, 
has the two types of transitions u —>■ u': 

(i) transitions inside the A/'-shell (close transitions), for which energy differences are 
Puu' = (ii ~ epA''^^^, and the maximal value L of the matrix element i^^^, is related to a 
transition between states of a j-shell; 

(ii) transitions between shells with major quantum numbers Af and A/'± 2 (distant tran- 
sitions) with Pj,,,/ = (i2 ~ EpA'^/^ and i^^^, ~ LA'^^^. 
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For the state 1, there are the three groups of terms in the sum (|3.11D, which are classified 



according to different combinations of the close and distant transitions in the product of the 
four matrix elements Those involving four close transitions have all the dimensionless 
differences x^^t ~ 1 and consequently F ~ 1.0 Thus, the contribution of these terms to 
the sum ( 3.11 ) is of the order L^. For terms with four distant transitions {x^^i ~ 



have F ~ A^^^. However, this large factor is compensated by the product of small matrix 
elements The same compensation takes place in the remaining terms with two close 
and two distant transitions, for which F ~ A^/'^. Therefore, the contributions of all terms 
in the sum ( p.ll|) are of the same order of magnitude L^, and we can make the following 



estimation: 

E^?2^23^3/41'5(^1-^f) = ^ Y.{i%,6{e,- S J.) . (3.19) 

Calculation of the last sum within the Tomas-Fermi approximation gives 

(^s-^J n{r)pUr){y' + z'fdv. (3.20) 



In this calculation we used the procedure described in Ref. |^T|, which includes averaging 



over the direction of the nucleon impulse and an utilization of the ansatz 

Y^^\(v)^,{r)5{e,-e,) = 4--n(r), (3.21) 
1 Pfk''^) 



where ?7,(r) = Cp^{r) {C = const) is the nucleon density and Pf{^) = y 2M[ej? — [/(r)]. 
Comparing ( |3.2(]| ) with the rigid-body moment of inertia 

^rig = / n{r){y^ + z^)dr, (3.22) 

we obtain Ps ~ '^rigippR/ ~ '^rigi^jp/^)'^, where R is the mean square radius of a 
nucleus, pp and jp ~ A^^^ are the mean impulse and the mean angular momentum of a 
nucleon on the Fermi surface. Thus, the parameter Ps has the order of magnitude h'^A^^^^/e^p. 
With these results we can get from ( |1.4] ) the perturbation parameter, jSgUj'^/a ~ {fiujjF/ 
An application of the perturbation theory implies that this value is small, i. e., the Coriolis 
interaction is smaller than a two-quasiparticle excitation energy. One can say that the 
perturbation theory is valid for adiabatic rotation. It is clear from Eq. ( |1.7| ) that Bs ~ SpA'^ 
and Bs/ A ~ A~^f'^. The above estimations refer to ND nuclei in the ground state where 
pairing correlations are reasonably strong, A ~ epA'"-^/^. 



4 The model of anisotropic oscillator potential 



In order to obtain quantitative results, we model the real self-consistent nuclear field as the 
axially deformed oscillator potential with the frequencies uj^ along the symmetry axis and 
(jjx in the perpendicular plane: 



^The necessary estimation for the amplitudes of the mriform pairing field, Di ^ D2 



can be obtained from Eqs. (A. 4) and (A. 



(4.1) 

[In(d2/A)]"^^l, 
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The use of this simphfied model is justified by a possibihty of deriving an exact analytical 
expression for the parameter (3^- It is known also that the model reproduces the experimental 
values of the rotational constants A and B for ND bands. Therefore we can expect the model 
to provide an useful insight in the rotational regime at superdeformation. 

In an anisotropic oscillator potential the matrix element is non-zero for four transi- 
tions. The two close transitions have the energy differences di = ±h{ujx — uj^), whereas the 
distant ones have d2 = ^h^u^ + tu^). The corresponding dimensionless parameters are 

where hcuo = 41A~^/'^ MeV. Hereafter we use the axis or frequency ratio k = b/a = cux/^z 
and the volume conservation condition uj'^uJz = "^o- It is also convenient to substitute the 
operator with its time derivative 

L = M{u;l-u;l)yz, (4.3) 

which has the same selection rules. 

For the fixed state 1, the sum ( p.ll[ ) involves 6 terms with the four close transitions, 6 



terms with the four distant transitions, and 24 terms with the two close and two distant 
transitions. The products of four matrix elements are equal with the semiclassical accuracy 
for all 36 terms of the sum, ^12^23^34^41 ~ ^(^z)ii- Summation of all these terms gives 

= 18(2A)Vf E(4)ii5(^i-^i^)- (4.4) 

The function is the sum of the functions F corresponding to all 36 combinations of the 
close and distant transitions. The sum can be simplified by using the symmetry properties 
( p. 18 ). It is convenient to represent the resulting function in the form 



$n(^^i, jy2) = ^(^^1, jy2) + 2Dln{vi, 1/2), (4.5) 

where 

V2) = fi^i, -^1, ^1, -^u 0, 0) + 2/(z/i, z/i, -z/i, -z/i; 2z/i, 0) 
+2(z/i/z/2)^[2/(z/i, z/2, -z/2, -ui; Vi + z/2,0) + 2/(z/i, -z/2, z/2, -I'l] I'l - «^2,0) 

+ f{^l, ^2, -^1, -^2] ^1 + ^2, -^1 + ^2) + f{Vl, -1^2, -Z^l, 1^2; - 1^2, -t^l - 1^2)] 

+ {l^l/l^2Y[f{^2. -V2. V2. -V2] 0, 0) + 2/(z/2, Z/2, -Z/2, -Z/2; 2z/2, 0)], (4.6) 

and the amplitudes Dx, D2, and the functions 7i are determined by Eqs. ( [A.5| ), ( [A.ll| ), and 
( A. 12 ), respectively. The two terms in ( [4.5|) describe the two distinct effects of the Coriolis 



force: the rotation-quasiparticle interaction and the modification of pairing. 
In the Tomas-Fermi approximation, we have according to the ansatz (|3.21|) 

Y,{t%,5{ex-eF) = SCM [ tl{v)^2M[eF - U,,,{v)]dv. (4.7) 
1 •' 

Combining the result of integration with the expressions for the rigid-body moment of inertia 
( pj.22D and the mean level density near the Fermi surface 

PF = j dvY,^\{v)^x{r)6{ei-eF) = SCM J pM[eF - f/osc(r)]cir, (4.8) 
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we obtain 

In the case of a normal deformation, the quantity k is close to unity and consequently 
V2'^ v\- Inserting ( [4.9|) into ( [4.4|) and using the above approximations we get for ND bands 



^ri g'^nc(t^l, ^2 



/9.(ND) = ^:ii_!ll_L^. (4.10) 



The function $„c approximates in "the close transition limit": 

$nc(^^i, ^^2) = J^ci^i, 1^2) + 2Dlnciiy2), (4.11) 

where 

^c(i^i, 1^2) = f {1^1, -1^1, 1^1, -I'l] 0, 0) + 2/(z/i, z/i, -z/i, -z/i; 2z/i, 0). (4.12) 
$n depends on z/2 only through the amplitudes Di, D2, and the function 

7^,(z/2) =8 + 81n2z/2 + ln4z/2. (4.13) 



The level density near the Fermi surface can be obtained by combining ([4.8|) with the number 
of nucleons A = J n{r)dr. The result, pp = SA/ep, depends on the Fermi energy which is 
found from the volume conservation condition a'^b = R^. The expression (|4.1CI|) has been 



obtained in Ref. ITT 



Figure 2a shows that $„c approximates the exact function very well. It shows also 
that the contribution of the rotation-quasiparticle interaction is small compared to that of 
the pairing modification. This result is explained by the interference of the two effects: 
the Coriolis force and the nonuniform pairing field A^^^. Neglecting the latter results in 
compatible contributions of the two terms of Eq. ( [4.11| ) as it is seen from Fig. 2b. This 
result is consistent with the Marshalek calculations 



For superdeformed nuclei the parameters ui and z/2 are both large. Thus, we should 
expect decrease of f3s. It is convenient to rewrite Eq. ([4.4]) by introducing, according to 
( [4.2| ), the new parameters ^ and k instead of ui and 1^2 . We express the sum ( [4.7| ) in terms 



of the rigid-body moment of inertia and the number of nucleons A: 

>^(4)ll5(£l-£i.) = 5(1 + p)3^2 ^O^rig- (4-14) 



We have now: 

/5^(SD) = . '^2.^, \ Z.. ^riM^ k), em k) = $„(//!, z/2). (4.15) 



e/'ji + ky ^3 

15/l'(l + P)3A2^"*5" 



The function $, along with its limiting case of the uniform pairing, is shown in Fig. 3. It 
can be seen that nonuniform pairing reduces /3s (SD) even more than /^^(ND). 

Taking k = 2 and A = 0.5 MeV (^ = 0.065) as the representative parameters for SD 
bands, we find from Fig. 3 that $ ~ 1. This yields the following estimation: /5s (SD) ~ 
h'^iA/epf ~ /?s(ND)A-2/3 and ~ A'^. The last estimation is correct for the SD bands 
in the A ~190 mass region where B/A ~ 10~^. Yet it overestimates the experimental value 
for 8^Zr(l) (~ 10-5) and i^^Gd(l), ^^^I)y{1) (~ 10"^). The later two bands have the smallest 
value of this ratio among all SD mass regions. Thus, a high deformation and nonuniform 
pairing do not solve the problem of the SD band regularity. 
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4.1 Limiting cases 

The limiting cases most interesting to us are: strong pairing, uncorrelated nucleons, and 
extremely large deformations. Reference [|lT| considers the limit of small deformations. 

For a very strong pairing (A S> fiuo), the size of the Cooper pair Rhuo/A becomes much 
less than the nuclear radius R. Rotation of such nucleus is described by the hydrodynamic 
equations of the ideal liquid P], according to which the second inertial parameter vanishes. 
In the following analysis, the key aspect is the nonuniform pairing. For strong pairing, the 
quantities A(i) and A(2) are proportional to A because, as follows from Eqs. (|A.4| ) and 



Di ~ A^ and D2 ~ A^. Therefore this limit is instructive since it allows to check the 
solution of the integral equation for A*^^-'.^ 

In the strong pairing limit the parameters i>i and 1^2 are small. It is possible to simplify 
the function F by expanding g{x^i,f) and h{x^^') in power series of x^^f, and then approximate 
with a necessary accuracy by 

F = V2 + D^Vi + DlVi + DlVe + D^Pg, (4.16) 

where P„ is the polynomial of the nth power in Xi,i,>. With this function, performing a 
calculation similar to the one we used to obtain we find the limiting value 

lim$>„.,) = -?r^^^^y. (4.17) 



A^oo 3 \UJj. + UJ 

Combining this result with Eqs. ( [4.4|) and ( |4.14| ) gives 



Thus, the parameter (3s vanishes in the hydrodynamic limit. 

The rotation of a very elongated nucleus with k = b/a ^ 1 exhibits some interesting 
physics. For this limit, the parameters and 1^2 are approximately equal 

ui = u ± Su, V = 5v = (4.19) 

^ 2A 2 A ^ ^ 

where z/ 3> 1 and Su/u = a/b ^ 1. The nonuniform pairing is also important in this case 
because the small amplitude ^ 1/z/^ is compensated by the large value of z/^. As a result, 
the kinematic moment of inertia, which is the sum of the standard cranking-model term and 
the Migdal one, is close to the rigid body value: 



Ck(l) _ CTi: . 



l-(^VHln2. 



Z/2 



(4.20) 



For the first approximation ui = z/2, the function ( [4.5|) vanishes, $„(//, i/) = 0. The next 
term of its expansion in gives the estimation of the second inertial parameter jSs ~ [a/by. 
We can say that a strongly elongated nucleus in the superfluid phase has the rotational 
regime which is close to the rigid-body rotation. The deviation are of the order (a/6)^. 
The physical interpretation of this phenomenon is straightforward: all nucleons of a needle 



^The solution for A'-'^^ has been verified in Ref. [ pf| by obtaining the hydrodynamic moment of inertia. 
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shaped nucleus with exclusion of the small sphere of the radius a in its center are involved 
in rotational motion. 



Let us now consider the normal phase. The right-hand side of Eq. (|3.11|) vanishes in 
the limiting case A = 0. This result is an artifact of the semiclassical approximation used 
in deriving the expression ( p.ll ). The correct formula obtained from Eq. (|3.4| ) with the 



limiting values of the Bogolubov amplitudes {u^ = 0, v,^ = 1 for p,y = I and Uy = 1, Vy = Q 
for = 0, where is the nucleon occupation numbers) is 



(sp) 



This expression describes the effect of the Coriolis force on single-particle motion. It will 
be shown in the next section that the cancellation of the leading terms in the sum of ( [4.21|) 
substantially reduces this quantity compared to 



4.2 The second inertial parameter for uncorrelated nucleons 

In this subsection we estimate the parameter (3 in the normal phase. In view of the cancela- 
tion mentioned above, we have to take into account the centrifugal stretching effect, which 
happens to be the same order of magnitude as ( [4.21| ). As with the superfluid phase, we 
will use the Green's function technique.^ Our starting point is the equations in the Hartree 
approximation]^ 

[e - Kiv) - V(r)]G'(r, r', e) = 5{v - r'), (4.22) 



where /i^ is the cranked single-particle Hamiltonian ( p.l| ) with the oscillator potential ( [4.1D 
and 

V(r) = j dciv,{v, q) j ^G(q, q, e) (4.23) 

c 

is the self-consistent potential. We assume that the two-body residual interaction Vq is the 
effective quadrupole one, 

^.(r,ri) = -f E(-l)"92M(r")g2-,(r;'), (4.24) 

where the quadrupole moment ^2^^ is defined in terms of the doubly stretched coordinates 

r = 1 — X +}—y + k — z. (4-25) 

UJq UJq UJq 

The interaction strength x is determined in a self-consistent way as follows: 



This interaction provides the full self-consistency for deformed nuclei [2S 



^There is an alternative method based on minimization of the system energy in the rotating frame as 
a function of the osciUator frequencies and the rotational frequency lo under the constraint of 

constant volume. For the fixed occupation of single particle states, this method gives the same result as the 
one obtained below. 

^This approximation is adequate for a separable two-body interaction we will use. 
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As usual, we proceed to treat the cranking term V with the perturbation theory by 
expanding the Green's function and the self-consistent potential in the series 

G = Go + Gi + G2 + G3 + ..., V = V(°) + V(^)+V(2)+V(=^) + ... . (4.27) 

The unperturbed Green's function is 

Go(r, v„e) = Y: G,(e)^,(r)(p:(ri), = ^ _ ^ ^^^^ _ ^^^y (4.28) 

with 6 —>■ +0. The occupation numbers refer to a nonrotating nucleus. We may notice 
that under the self-consistent condition 

^x'^x = ^^x^y = ^z^z, ^x,y,z = '^{''T'x,y,z + ^/'^)uPu (4.29) 

V 

{nx,ny,nz are the oscillator quantum numbers) V'^"-' = 0. Thus, the average potential is 
modified by rotation only. 

The odd corrections to the self-consistent potential are lacking, V*-^*^^-* = 0, due to the 
different symmetry properties of the operators q2fi and ix under the time reversal. Conse- 
quently the third order correction to the Green's function is expressed as 

= GoVGqVGoVGq + GqVGoV'^^^Go + GqV^^'^GoVGo. (4.30) 

The first term generates the interaction of rotation with single-particle motion. It yields the 
quantity P^^^ ( [4.21| ). The last two are responsible for the centrifugal-stretching effect which 
is described by the expression 



o r rip 

/^r^ = --,f ^tT{V^'^Go{e)VGo{e)VGo{e)}, (4.31) 



V('Hr)=-xE(-l)''^2-M(r)/;^tr{g2M[Go(^)VGo(5)VGo(e) + Go(£)V(2)Go(e)]}. (4.32) 



j de 

c 

where the correction to the mean potential V^^-* is obtained from the equation 
The solution of this equation has the form 



V(^)(r) = -c.^ET-T;;-(-l)"Q2>2-,(r), (4.33) 

fj. J- + ^ 

where 

= E 1(92^)12^^^^, Q£ = / ^tT{q2,Go{e%Goie%Goie)}. (4.34) 



The last quantity is the second correction to the nuclear quadrupole moment due to rotation. 
Its explicit form is 

Qg = E(g2.)l2^23^3lEA| . , (4.35) 
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where the permutation of indices i/ = 1, 2, 3 by the operator Pi is subject to the rule i/mod3 = 
u. It is obvious that the non-zero corrections have the components with /z = 0, ±2. The 
denominator in the sum (|4.33|) renormalizes the interaction strength. The straightforward 
calculation of and the use of Eq. ( [4.26| ) with the zero-order density matrix p gives 



+ X'^m) ~ 2x- Combining ([4.31|) with (|4.33| ), we have 



= S (4.36) 



We can now calculate the two contributions to the parameter /?„ by summing over the 
quantum numbers n^, Uy, and n^. The anisotropic oscillator potential allows to find an exact 
solution of the problem. At first we find the corrections to the quadrupole moments 



(2) _ 5 hj:,{i-5k') 



Then by using (|4.36| ) we obtain the contribution of the centrifugal stretching effect 



= 3^3fc4(p_i)2 (^^ - + 76fc^ - 15^ + 1). (4.38) 



Finally, after some fairly tedious calculations of the sum (|4.21|) we get 



= ^^,J^;_^y ik' - lOk' - Uk' - lOk' + 1). (4.39) 
Adding the last two quantities gives us the parameter (3 in the normal phase: 

5%^k^-ioe + i 

^" 6u;i A;2/3(A;2 + 1) ' ^^■^^> 
if we use the following formula for the rigid body moment of inertia: 

ujzk"^ ' 

Parameter /5„ is substantially reduced compared to /3s, /5„ ~ /3s(SD)A^^/^. This 

can be explained by canceling main terms in the sums (|4.21| ) and ( [4.35| ). That is exactly why 
the corresponding values (3^^^ and Q^2ii proportional to S^. Such result is predictable 
because the Hamiltonian for the anisotropic harmonic oscillator can be diagonalized 
exactly pOl- Its eigenstates are characterized by the number of rotating bosons. To find 



^rig = -T|(fc' + l). (4.41) 



(3^^^ and we have to calculate first the expectation values of the operators (.r^ and in 
this rotating basis. Then these quantities must be expanded in powers of uj. Because these 
operators are represented by quadratic forms in the rotating bosons, their mean values and 
therefore all the terms of the series are proportional to the linear combination of S^,., S^, 
and S^. 

Another peculiarity of the solution (|4.40|) is that /3„ < for the prolate nuclei with 



1 <h/a< 3.15, whereas (3s is always positive. The formal cause of this effect is a negative 
value of (3^^^^ and the inequality \(3^'^'^\ > (3^^'^'^ >0 which is fulfilled for the above indicated 
deformations. In the superfluid phase, the term responsible for the rotation-quasiparticle in- 
teraction may be also negative, but it never exceeds the contribution of the pair modification 
effect (see Fig. 2a). 
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5 Analysis of experimental data 



We have shown in the preceding section that the second inertial parameter B is negative 
in the superfluid phase and positive in the normal one. The two hmiting cases allow us to 
reconstruct the B{I) dependence for the parametrization ( |1 . 1| ) of a rotational sequence with 
(vra) = (+0). Comparing the formulas ( [4.10 ) or ([4.15|) with ( 4.4C ) we conclude that the 
ratio B/A has to change sign with increasing the spin / in a band and approach its limiting 
value Bn/An ~ A'^/^ for high I. 

The limiting ratio for a real nucleus can be obtained from Eqs. (|4.4CI| ) and ( [4.41|) if we 
suppose that the r.m.s radius and the deformation are exactly the same for neutron (z/) and 
proton (tt) systems. The first condition implies that oscillator frequencies of neutron and 
proton potentials satisfy the relation c^or = ^^oC^At/A) (r = tt, z/ and A^ is the number of 
nucleons of a given type). The second one results in the identical ratio of the frequencies 
along the principal axes for the both potentials: 



UJyr '■ Ulz 



m : m : I. 



(5.1) 



For integers m and /, the states with the same number of quanta Afmi = m{nx + riy) + luz 
form a deformed shell. Assuming that, for a given number of nucleons A^, all the A/'mrshells 
are filled, one can express the sum Tj^t in the form 



^ZT^ZT — ^Ori'^XT^yT'^ZT 



,1/3 



1/3 



(5.2) 



The above formulas allow us to derive the ratio B/A for a nucleus consisting of Z protons 
and neutrons in the normal state: 



Bn ik^ 
-^ = -3.205^ 



10fc^ + l)p/3 

(F + 1)3^8/3 



a) + 7 



(5.3) 



This result holds for a nucleus with an arbitrary deformation k = b/a = m/l. 

We concentrate first on SD bands. Most of them are not connected to lower-lying states of 
known excitation energy, spin, and parity. Thus, their exit spins Iq are unknown. Tentative 
spin assignment is used to take advantage of the formulas ( |1.2| ) to find the experimental ratio 
B/A. To analyze this quantity, we will take into account two basic ingredients: shell gaps, 
which stabilize the shape, and intruder orbitals involved in the alignment. The nucleon- 
configuration assignment of a band is generally based only on the behavior of the dynamic 
moment of inertia and the quadrupole moment in a given band. The last quantity. 



Qo = 6.05 ■ IQ-^A^/^ 



A;2/3 



-eh. 



(5.4) 



remains remarkably constant as a function of spin within a band. This proves that the 
deformation k remains practically unchanged as / increases. We use the experimentally 
observed value of Qq to find the axis ratio 6/a, which is required for calculation of the 
quantity ( ^^ . 

The ratio B/A extracted from the measured energy of 7-transitions in the four SD bands 
of the A = 150 mass region are shown in Fig. 4. The parity and the signature of these 
bands is assumed to be (+,0). We also use the adopted spins for their lowest levels. The 
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band ^^^Dy(l) belongs to the doubly-magic nucleus with the proton Z = 66 and the neutron 



N = 86 gaps in single-particle spectra at the same deformation |31|. The gaps decrease a 



level density and considerably reduce the neutron and proton pairing correlations. There 
is no direct experimental indication of pairing correlations in this band. The theoretical 
calculations |32| show that their inclusion leads to a better description of the kinematic and 
dynamic moments of inertia, pairing correlations being more important at the low spin range. 
The plot shows that there are two distinct regions in the variation of the ratio B/A versus 
/. The lower part of the band exhibits a sharp increase of this quantity. It then changes 
sign at the spin Ic = 36 and approaches the plateau value of (|5.3| ) at the top of the band. 
Such behavior of the ratio apparently shows that the static pairing correlations of neutrons 
and protons are quenched simultaneously. This fact also proves that most part of the band 
belongs to the normal phase. 

The band ^^'^Gd(l) is one of the few examples of SD bands which exhibits backbending. 
The 7r6^ pair alignment opens up the proton shell gap Z = 64 at the same deformation as 
the neutron shell gap = 80. Thus, above the backbending this band becomes similar to 
the doubly- magic ^^^Dy(l) except the gap Z = 64 is less pronounced than the one at Z = 66. 
Besides, the neutron gap = 80 is somewhat smaller than = 86. These factors enhance 
a level density and favor pairing correlations. As seen in Fig. 4 the behavior of B/A for this 
band in the low-J region is the same as that for ^^^Dy(l) if we scale the axis of ordinates by 
the factor two. Accordingly, the critical value is somewhat larger, Jc = 38. 

The features observed at low spins in the dynamic moment of inertia of the band ^^°Gd(l) 
have been explained in terms of consecutive alignments of the and 7r6^ pairs |^|. For 



the configuration 7r6^z/7^ all levels below the Z = 64 and A^ = 86 shell gaps are occupied. 
The former is found at slightly smaller deformation than the latter. This factor diminishes 
the neutron gap and enhances neutron pairing. The Woods-Saxon and the relativistic 
mean field |^B[ calculations make evident that static proton and neutron pairing exist at low 
spins, / < 48 {huj < 0.55 MeV). In addition, the full self-consistent HFB calculations with 
the particle number projection show that the effect of pairing on the moments of inertia 
in ^^''Gd(l) is about twice as important as in ^^^Dy(l). It is apparent from Fig. 4 that the 
static pairing correlations in the band ^^°Gd(l) are stronger than in ^'^^Gd(l). One should 
expect even more stronger pairing correlations in the newly discovered prolate deformed 
band ^^^Er(2) |^4[, since the proton Fermi level at Z = 68 lies in the region of a high level 
density above the Z = 66 gap. The experimental B/A dependence of Fig. 4 is consistent 
with this prediction. It is seen that this ratio does not exceed the value —5 x 10~^ and 
does not show the plateau.j] The bump seen at / = 44 (cu = 0.57 MeV) can be attributed 
to the alignment of a pair of ^13/2 protons in agreement with the calculations of Ref. 
Thus, the plots of Fig. 4 show the correlation of the spin dependence of the ratio B/A with 
the level density near the Fermi surface: the higher the level density, the stronger pairing 
correlations, and the less marked the plateau. 

The high deformed (HD) bands in the A = 190 mass region are related to the Z = 80 
and A^ = 112 shell gaps. Most of these bands have similar values of 3*^^^ that exhibits a 
smooth rise as a function of rotational frequency. This rise is attributed to the gradual 
alignment, in the presence of static pairing correlations, of 213/2 protons and ji5/2 neutrons. 
The calculations with pairing are able to reproduce the general trend seen in experiment. 
The bands ^^^Hg(l) and ^^''Pb(l) are of central importance because their spins, parities, and 

''it is worth mentioning that the nonaxial band ^^^Er(l) demonstrates the B/A dependence with the 
critical spin Ic = 31 and a long plateau. 
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excitation energies are known P5| BBl. The plots of Fig. 5 for these bands demonstrate the 



gradual rise of the B/ A ratio that confirms the presence of static pairing correlations. 

Now we turn our attention to ND bands. There are several bands of isotopes Er, Yb, 
and Hf in which the static neutron pairing gaps are predicted to collapse. However, the 
proton system still has strong pairing correlations. Accordingly, the B/A vs. I plot for 
these bands exhibits a sharp rise, but does not approach the plateau. The yrast band of 
^'^Zr is an exception. Because protons and neutrons occupy in this nucleus similar orbitals 
near the Fermi surface, quasiparticle alignments and the elimination of pairing gaps occur at 
similar spin. Besides, the deformed shell gaps at Z = = 38 and the low moment of inertia 
favor the transition to the normal phase. Combination of these factors makes the pairing 
properties of the ND band of ^^Zr similar to those of SD bands in the A = 150 mass region. 
Fig. 6 plots, as a function of spin J, the B/A ratios for the ND and SD bands of this nucleus 
determined from the data of Refs. and [^, respectively. In the SD band ^'^Zr(l), static 
pairing correlations are quenched completely due to the high rotational frequency. This 
inference is supported by the coincidence of experimental points with the plateau Bn/An- 
It also becomes apparent from this figure that at high spins the experimental ratio for the 
ND band reaches the same plateau. The difference in the limiting value of Bn/An due to 
the difference in deformations (/5 = 0.43 and 0.55 for the ND and SD bands, respectively) is 
insignificant. The low spin part (/ < 18) of the ND band is compatible with the transitional 
nature of the 7-soft nucleus: small f3 and a noticeable triaxiality. The alignment of the two 
5^9/2 quasiprotons and the subsequent alignment of the two gg/2 quasineutrons are clearly 
seen in Fig. 6 as the humps A and B. Beyond the second alignment, a striking change in 
deformation occurs in the interval of spins J = 18 — 22. After the spin J = 24 the rotational 
behavior is compatible with the rigid rotation of a high deformed axially symmetric nucleus 
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The characteristic behavior of the ratio B/A with the critical spin 1^ and the pronounced 
plateau have also been found in the SD bands of the A = 150 nuclei having configurations 
different from (+,0). The bands ^^^Th{2) [the 7r[301]l/2 hole in the ^^^Dy SD core] and 
^^^IIo(3) [the ^^^Dy SD core coupled to the 7r[523]7/2 orbital] show the B/A dependence 
similar to that of the ^^^Dy(l) band. The pair of identical bands ^^°Gd(2) and ^^^Tb(l) have 
the dependence similar to one of ^^°Gd(l). All these bands have somewhat higher values 
of B/A than do their (+,0) counterparts. This proves that an odd nucleon or a particle- 
hole excitation reduces pairing correlations due to the blocking effect. Such phenomenon is 



characteristic of the static pairing regime ||l9[. The SD bands from different mass regions, 
^^^Ce(l), ^^^Ce(l), and ^°Zn(l), exhibit the same behavior of B/A. A strong configuration 
dependent effect is observed in the bands where the odd neutron is placed in the ji5/2 intruder 
orbitals. For such bands, the B/A ratio is positive for all spins. The examples include the 
bands i^9Gd(l) (configuration TreVT^), ^^^Dy(l) {7rQ^u7^), and ^^^Dy(l) (vre^z/T^). A single- 
particle degree of freedom seems to destroy the typical behavior of the ratio B/A. More 
efforts are needed to explain this interesting feature. 

These numerous examples prove the universality of the transition from the superfluid to 
the normal phase for SD and ND bands. This universality can be represented, according to 
Ref. [^, by the effective rotational Hamiltonian, 

Hes = al^ + (///c - l)bl^ + cl^ (5.5) 

which describes the states of the (+, 0) band in the transition region. The parameters a, 
b, c, and the critical spin are the subjects of the microscopic theory, which has to take 
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into account static, dynamic, and uniform pairing. Incorporating the critical spin, which 
can be found from the experimental plot B/ Ays. /, this concept of the superfluid-to- normal 
transition is free from ambiguities characteristic of the approach based on a change of the 
single-particle spectra 0, ||. 

Using the results of our analysis we can now explain why some SD bands have extremely 
regular rotational spectra. Figure 4 shows that the most part of the bands ^^^Dy(l) and 
^^^Gd(l) belongs to the plateau with the ratio of the inertial parameters B/ A ~ 10~^. So 
does the whole of the SD band ®^Zr(l), for which this ratio is 10~^. The plateau is the 
manifestation of the normal phase with the anomalous small ratio (|5.3| ), Bn/ An ~ A~^^^ . 
The above values for the SD bands agree with this estimation. Therefore the extreme 
regularity is explained by quenching the static pairing correlations in the lower parts of 
these bands. On the contrary, the transition in the yrast band of ^^Zr occurs in its upper 
part. Accordingly, the top of the band has the same properties. It is important to note that 
the bands in which proton and neutron pairing gaps are present [^^^Er(2) and all the SD 
bands in the A = 190 mass region] and the bands with proton pairing alone [^^'^Yb(yr) and 
-'^^^Hf(yr)] are regular to a lesser extent. 

We should mention one more feature which requires further investigation. The down- 
sloping of the B/A dependence is observed at the top of ^^^Dy(l), ^''Zr(l), and other SD 
bands with extremely high spins. Because the quantity Bn/An is a decreasing function of 
the deformation k, it would be natural to explain this feature by the increase of the nuclear 
elongation due to the enormous centrifugal stretching at the end of these bands. 



6 Conclusion 

Despite the vast amount of data collected and various theoretical interpretations suggested, 
a detailed understanding of many properties of SD bands has yet to be achieved. Pairing 
correlations are just one example of such properties. The presence of static pairing in SD 
bands is usually established by studying the behavior of the dynamic moment of inertia Q"*^^-' 
as a function of the rotational frequency u. A band crossing associated with a quasiparticle 
alignment leads to a impressive decrease in ^J*^^^ with or a hump in this dependence. This 
gives an indication that static pairing correlations are present in that part of a band where 
such irregularities occur. 

In this paper, the investigation of pairing correlations is based on the spin dependence 
of the second inertial parameter B. This quantity, which is proportional to the difference 

CJ(I) _ 

the high- J limit, turned out to be a more sensitive measure of the change 
in pairing correlations than ^^"^K The new method requires spin-signature assignments of 
the band states. However, it gives more definite information about the superfluid-to-normal 
transition in a band. The most important results obtained in this paper can be summarized 
as follows: 

(i) The exact semiclassical expression for the second inertial parameter in the superfluid 
phase has been found by taking into account the effect of rotation on the Cooper pairs in the 
gauge-invariant form. The presence of nonuniform pairing reduces the nonadiabatic effect 
of rotation. Its influence increases strongly at superdeformation. The nonuniform pairing 
allows one to find correctly the interesting physical limits for the second inertial parameter. 

(ii) The limit of zero static pairing is of special interest. It permits the function B{I) to be 
reconstructed by interpolating between the values of B in the superfluid and normal phases. 
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The anisotropic oscillator model calculations show that there are two distinct regions in the 
variation of the ratio B/ A with I. The lower part of a band is characterized by a gradual 
decrease of pairing. Accordingly, being negative the ratio B/ A exhibits a sharp increase. It 
then changes the sign at the spin and approaches the positive value characteristic of the 
normal phase. The critical point Jc, B{I(^ = 0, is a signature of the superfluid-to-normal 
transition. The transition manifests itself in the modification of the rotational spectrum of 
a band. 

(iii) The experimental spin dependence oi B/ A indicates an agreement with the theoret- 
ical prediction and demonstrates the universality of the transition to the normal state. This 
agreement is not a trivial fact because our calculations are based on the simplest model of 
a nuclear potential and do not take into account pairing fluctuations in the normal phase. 
Nevertheless, the agreement is not accidental because the universal dependence of B/ A on 
/ has been observed for a large number of SD bands and some ND ones. 

(iv) The universal spin dependence oi B/ A explains the extreme regularity of some SD 
bands. The characteristic feature of this dependence is the pronounced plateau in the upper 
part of a SD band (/ > /c) corresponding to the normal phase. The calculated ratio in this 
part of a band is extremely small, B/ A ^ A^^^'^. Thus, the closer the critical point Ic to 
the exit spin Jq, the more regular its rotational spectrum. The spectacular examples are the 
bands i^4Gd(l) and ^^2Dy(l) having B/A ~ 10"^ 

(v) Some new features have been observed in the upper parts of SD bands. The investi- 
gation of this region, which is free from pairing correlations, is extremely important for our 
understanding of the microscopic structure at the superdeformed minimum. 
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Appendix A. Solution of the integral equations for the 

nonuniform pairing field 

The effect of rotation on pairing correlations is described by the first, A'^^)(r), and the second, 
A^^^(r), corrections to the pairing field, which enter into Eq. ( ^.41) . We have seen in Section 
1^ that the integral equations that have to be solved are of the general form (2.30) and ( ^.31 ) 



for even and odd corrections, respectively. It is convenient to introduce into these equations 
the operator Cx that is a function of the space coordinates only. 

We start by considering the equation for A*^^). Using the relation = i{pi —^2)^12 we 
get from (|2.31|) for i = 



1,2 



The equation is satisfied if we assume that 
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where the amphtude Di is determined after substituting ( |A.2| ) into Eq. ( |A.1| ) and integrating 
over r: 



■X 1 2 
12 I 



1,2 



dp 



1 



27r (p2 + A2)(p2 + A2 



:c 1 2 
12 1 



1,2 



(Pl -P2)^ 



C" 



27r (p? + A2)(p2 + A2) 



(A.3) 



The solution transforms to 



Si, 2 ^12 


^5'(a;i2)5(£i-£F) 


Z)l,2 Kl2 


^a;i25'(a;i2)5(£i-£F) 



(A.4) 



in the semiclassical approximation. The summation over the state 1 is to be understood 
as the integration over its quantum numbers. For the anisotropic oscillator potential, the 
amplitude can be expressed in a simple analytical form 



^1 



(A.5) 



The function g{x) and the parameters vi and V2 are determined by Eqs. ( p.l6| ) and 
respectively. 

The equation for A*^^)(r) after introducing f^. becomes 



E 

1,2,3 



dp 



Q3{P,Pl,P2,P3) 



C 



12 "23 



2n (pf + A2)(p2 + A2)(pi + A2 

dp {pi - + 4A^ 
; 2^(pf + A2)(p2 + A2) 



0, 



(A.6) 



where the polynomial function Q3 of the third order in {p,Pu) depends also on the amplitude 
Di. We try to solve this equation by making substitution 



A(2)( 



2, ,2 



4A3 



D2t 



(A.7) 



Applying the same procedure as before, one can find the amplitude D2 in the semiclassical 
approximation 

El,2,3 ^?2^23(^x)310(a;i2, X23)S{ei-eF) 



j:i,2,3m)i2mx,2)S{e,-eF) 
where the function h{x) is determined by ( |3.16| ) and that of has the form 

1 



(A. 



■[—xy{l — Dix ){1 + xy — Diz )g{x) 



2x'^y'^z'^ 

+ y^{l - Dix' - Diz^)hiy) - yz{l - Diz^){l + yz- D,x^)g{z)]. 
Their symmetry properties are 

0(2:, y, x) = y, z), -y, -z) = y, z). 



(A.9) 
(A.IO) 



In the oscillator potential, the sum over the states 2 and 3 in the numerator of (|A.8|) comprises 
16 terms including four with two close transitions, four with two distant transitions, and eight 
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terms with one close and one distant transitions. Performing summation in the semiclassical 
approximation in the numerator and in the denominator of ( [A .81) , we find 

D2 = [40(z/i, Z/i - Z/2, -Z/2) + 40(z/i, Ui + U2, U2) 

+ 2i/i, ui) + 0(z/2, 2z/2, z/2) + 40(z/i, 0,-z/i) + 40(z/2, 0-iy2)]n'\ui, U2), (A.ll) 
^(^^1, 1^2) =8 + Ah{ui - z/2) + 4/i(z/i + z/2) + /i(2z/i) + h{2u2). (A.12) 



where 



For the monopole pairing interaction, the pairing field is uniform and the first correction 
A*^^^ vanishes. The coordinate independent solution for the correction A*^^^ can be found 
from Eq. ( |A.6| ) after averaging over r. The resulting expression may be written in terms of 
the kinematic moment of inertia, 



a;2 a^(i) 



2pF OA 



(A.13) 



1,2 



in agreement with the result obtained by Marshalek ||T3[. From the theoretical viewpoint, 
this solution is not correct because it violates the current conservation. 



Appendix B. Calculation of integrals 

In this Appendix we give a brief outline of the technique used in calculation of the integrals 
which are necessary for obtaining the function F (|3.12|) and for the solution of the integral 
equations ( |A.1| ) and ( |A.6[ ). All the relevant integrals can be done exactly with the method 
proposed by Feynman in the quantum electrodynamics pSf. The method is based on the 
identity 

1 = (n-l)\ [ dh [ dt2...7 , ^^^^ , r^, (B.l) 

aia2...a„ 7 J J [aitn-i + a2(in~2-in-i)+ ••• +an(l-ii)J" 

which is proved by a direct calculation. 

The simplest integral is the one involved in the sum ( |3.1U| ). It is solved by using ( |B.1D 
as follows: 

J ^ fdPi 1 

' J 27r (pf + A2)(p2 + A2) 

J 2n [pl + phQit)]^ 2pi2J Q{t) 2A2^12aJ' ^ ^ ^ 







where 



g(t) = -t' + 1 + (5^ (5 = A/pi2, Pi2=Pi-P2- 

Four integrals appear in the first sum of Eq. (|A.6|) . All those are of the same type. As an 
example, we consider the term proportional to the square of the amplitude Di. The relevant 
integral is 

dvi Q2ip;pi,P2,P3) 



J2 = A 



2n (pf + A2)(pi + A2)(pi + A2) 



27 



A [dt, Jdt, + + (B.3) 







where Q2 = P + - + P2P3 + ^ 

(5(^1,^2) = -[(1 - c')ti + 0% - If - (1 - c2)ti - 0% + l + 5\ 

and c = Pu/Pn- By making the substitution ui = {1 — c)ti + ct2, U2 = t2, we get after 
integration over dpi 



J, 



where Qi{ui) = —uf + (1 — c)ui + 6"^. Straightforward calculation of these integrals gives 

J2 = [Xugixu) + X23g{x23)]- (B.5) 

ZAX13 



Finally, let us consider the integral ( |3.9| ). It is solved in the same manner as the preceding 
ones. After using the identity (p.l| ) and integration over dpi, the substitution ti = [ui — 
{d — c)u2 — CM3]/ (1 — c?), ^2 = M2, ts = M3, where c = P12/P14, d = Pn/pu, leads to four triple 
integrals which may be solved without a problem. 
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FIGURE CAPTIONS 



Fig. 1. Relative deviation of energies E in the superdeformed band ^^''Pb(l) as compared 
with the ground state band of ^^^U and the ground vibrational state band of the H2 molecule. 
The deviation is calculated from the formula {E — Ej-ig) / E^g, where Ej-^g = AI{I+l), and the 
parameters A are found from the energies E^{4) of the 6^4 transitions. The experimental 



data are taken from Refs. and 



Fig. 2. Comparison of the functions relevant to the second inertial parameter for ND 
bands. Part (a) represents the function $„ (solid line), $„c (|4.11|) (dashed line), and 

their parts, the functions JF ([4.6|) (dotted line) and J-'c ( |4.12| ) (dash-dotted line) which describe 
the effect of rotation-quasiparticle interaction. Part (b) is for the limit of close transitions. 
The function $„c (solid line) and its constituents J-'c (dashed line) and 2D2TCC (dotted line) are 
shown for A*^^^ = 0. The axis of abscissas shows the dimensionless quantity ui corresponding 
to the close transitions, while that for the distant ones is fixed by the representative value 
i>2 = 10 for all the plots. 



Fig. 3. Plot of the function $, to which the second inertial parameter for SD bands ( [4.15| ) 



is proportional, against the dimensionless quantity ^ for the axis ratio b/a = 2. The solid and 
dashed lines correspond to the exact value and the limit of uniform pairing, respectively. The 
abscissa scale must be multiplied by a factor of approximately 7.7 for nuclei in the A ~ 150 
mass region in order to obtain the gap energy in MeV. 

Fig. 4. Ratio B/A versus spin for some SD bands of the A = 150 region. Expressions 
( |1.2|) are used to extract this ratio from experimental data taken from Refs. ^ and p4| . 
The error bars (if they are greater than symbol sizes) include only the uncertainties in 7-ray 
energies. The uncertainties in the spin assignment are immaterial for all bands [with the 
exception of ^^^Dy(l)], since the spin variation in 2h would merely shift the curves along 
the abscissa. The experimental points for the band ^^^Gd(l) are shown above the 7rii3/2 
backbending. The solid straight line is the ratio Bn/An ( |5.3| ) for the normal phase with the 
deformation b/a found from the quadrupole moment 



Fig. 5. The same as in Fig. 4 for the two SD bands of the A = 190 region. Experimental 
data are taken from Refs. [§, ^ ^ . 



Fig. 6. Ratio B/A versus spin for the yrast ND band (solid circles) and the SD band 



(open circles) of Zr. Experimental data are taken from Refs. ^Tj and ||38[. The solid 



straight line is the ratio Bn/An ( ^.3|) relevant to the deformation of the SD band. 
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